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We describe and test a method known in the literature as low mode averaging to improve Eu- 
clidean two-point functions in lattice QCD using the low-lying eigenmodes of the Wilson-Dirac 
operator D. The contribution from the low modes is averaged over all positions of the quark 
sources while the contribution from high modes is calculated in the traditional way using one 
source point per lattice. We apply this method to different baryon and meson two-point functions 
and we compare the improvements using the eigenmodes of the non-hermitian operator D and the 
eigenmodes of the hermitian operator Q = 75D. The convergence strongly depends on the parity 
of the states. 
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1. Introduction 



Solving the Dirac equation Dy/ = tj may be a very time consuming task at small quark masses 
as the condition number of D is proportional to the inverse of the quark mass. Some approaches 
involve the computation of the low-lying eigenvectors of 75D and use these to deflate the Dirac 
operator, see e.g. ref. [JTJ, §j. These eigenvectors can also be used to reduce the noise of the signal 
with a technique known in the literature ||] as low mode averaging (LMA). This consists of 
decomposing the quark propagator into a sum of high and low modes y = y/high + Vtow and in 
averaging the low mode contribution over all lattice points. We define the LMA two-point function 
as, 

Clma(0 = Q ow (t) + C pa (0 -00 , (1.1) 

where C pa is the exact point-to-all correlation function, calculated for a single source point. This 
definition satisfies (Clma(O) = (C pa (0)> however, the errorbars are reduced, due to more sampling 
per lattice. In the present work we apply the LMA technique to meson and baryon two-point 
functions with degenerate quark masses. We study the improvement as a test for future work on 
hadron spectroscopy and three-point functions. 



2. Mesons 

The present study is based on 100 configurations with lattice volume V = 16 3 x 32 generated 
with the quenched Wilson action at j3 = 6.0175 using Chroma [ph. This corresponds to a lattice 
spacing of a « 0.2093 a -1 / 2 « 0.093 fm. For each configuration we computed the lowest 30 eigen- 
vectors of the massive hermitian Dirac operator y$D\ui) = A,|w;) using the algorithm by Kalkreuter 
and Simma [|J at K = 0.1557 corresponding to m % 425 MeV and m n L « 3.2. We can recon- 
struct the contribution of the low modes to the quark propagator as [|7p, Dj^, = Y,i j-\ u i) ( u i\Y5 while 
Ciow(0 is given by, 

Ciow(0= E <K 7 (x)]r 5 r < Mj -(y) |y 5 r | M7 0)> , (2.1) 

where yn =X4 + t. 

Clearly the LMA technique works best when the low modes are dominant. This happens 
for the % as can be seen from fig. [I] (left) where we compare the smeared-smeared % point-to-all 
correlation function with the low mode contribution. However, for positive parity mesons (see fig. |] 
(right) for the ao), the low modes saturate the two-point function very slowly. This can also be seen 
by comparing the LMA and point-to-all smeared-smeared effective masses of fig. |2| for negative 
parity mesons with the positive parity ones of fig. [| In the case of the ao the low mode contribution 
even has a wrong curvature in the central region, see fig. [I] (right). 

The meson correlators differ in the T structure of the interpolating fields: for the 71 (the best 



case) we have a total T product y^F = t in eq. (2.1) while for the ao this reads y^F = 75. One 



way to get rid of a 75 factor is to calculate the right and left eigenvectors of the non-hermitian 
Dirac operator, D|r,-) = A,|r ; ), (£i\D = and to construct the quark propagator as [Q, Dj^, = 
E/ j: I r i)(^i I- m this way the meson two-point function becomes, 
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<7 pa (i) 




Figure 1: Low mode saturation (30 modes): n (J PC = + , left) and oq (J pc = ++ , right) two-point 
functions. 
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Figure 2: Effective masses: % (J PC = + , left) and p (/"- = 1 , right) 
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Figure 3: Effective masses: a\ (J PC = 1 ++ , left) and b\ (J^ = 1 + , right) 



jpc - 



ciow(0= I j-r^Wirir/WX^WiriryCy)). 



(2.2) 



Note that = 5,-y and it can easily be seen that due to the property = 75D75, A* is an 

eigenvalue whenever A,- is an eigenvalue, with a left eigenvector (£i\ = (r^ys and a right eigenvector 
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m 


"Jlma 




0.2003(49) 


0.2019(36) 


p 


0.3965(66) 


0.3952(41) 



Table 1: Fitted meson masses in lat 
tice units. 



We used the Arnoldi method implemented in the ARPACK 
library to compute the non-hermitian eigenvectors. This method 
is much slower than that of Kalkreuter and Simma. However, 
the eigenvectors are independent of the hopping parameter k 
and when changing the quark mass we only need to rescale the 
eigenvalues A. Unfortunately, we could not find any signifi- 
cant improvement using 15 pairs of non-hermitian low modes 
over the standard point-to-all method. This may be related to the fact that the ratios IA30I/IA1I 
were found to be by factors of approximately five smaller in the non-hermitian case than for the 
hermitian case. 

From the fit to the correlators we found a 30% improvement on the extrapolated % and p 
meson masses for the hermitian LMA, as displayed in table [l[ 
We also applied the method to the local axial- 
vector current (fig. |]), to test the improvement 
on the pion decay constant f % . Following the 



conventions of ref. ||10(] we set, 
m«fx = <0|7 4 7 5 |?r> = 



2k v Z- ■ (2.3) 



/App 



We denote the amplitudes of correlation 
functions with a smeared pseudoscalar sink and 
a local axial-vector source as A&p and with a 
local pseudoscalar source as A PP . The improve- 
ment of LMA is very striking in this case, 




Figure 4: Local axial-vector current. 
fx = 0.079(27) , fn.LMA = 0.082(12) . (2.4) 



3. Baryons 

The low mode averaging technique for baryons is more challenging due to the way the inter- 
polators are contracted. For the nucleon N, the A + and the A we use the following interpolators, 

N{x) = e abc u(x) a {u{x) T b Cy 5 d{x) c ) , (3.1) 
A + (x) = Babe [2 («(x)JCy M d(x)b) u(x) c + (u(x)lCYpu(x)b) d(x) c ] , (3.2) 
A(x) = e abc [2 (u{x) T a Cy 5 d{x) b ) s(x) c + {u(x) T a Cy 5 s{x) b ) d(x) c - {d{x) T a Cy 5 s{x) b ) u(x) c ] , (3.3) 

where C is the charge conjugation operator. Note that we only study the mass-degenerate case, 
m u = rrid = m s . In the following, we consider the contractions for the example of the nucleon 
two-point function, 

(N(y)N(x)} = {£ a bcea'b'ciCy 5 ) a p(Cy 5 ) al p l (P ± ) 7 Yd(y) b >i 3/ d(x) b pu(y)a'a'Hx) aa u(y) c yTl(x) cr ) , 

(3.4) 

where P± = i(l± 74) is the parity projector. 
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The expression eq. (3.4) can be written in terms of the left and right eigenvector components 
in a way that the source term is decoupled from the sink term: 



(N(y)N(x))=£ 



1 



{£a'b'd{Cys)a'p' [ri,b>,i3'(y) r jy,a>(y)rLc>,y(y)] } (P±) 



jk XiXjX k 

< {eabc(Cr 5 ) ap £* Kp (x) [£l a , a (x)t kiC>7 (x) -^ c , 7 (*)4, fl ,«to] } 



(3.5) 



For the hermitian case, y=,D\ui) = Xi\ui), = \uj) and = {uj\y=,. The cost of computing the 
above expression increases with the third power of the number of eigenmodes used. To optimize 
this computation, we split up each eigenvector r into twelve r a a spin-colour components and store 
them as lattice objects. Multiplications like r iaa rj b p are optimized in QDP++. With these and 
other optimizations the cost of the contractions for 30 eigenmodes becomes negligible relative to 
the cost of solving for the propagator. 
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Figure 5: Nucleon (J p = | ) two-point function and effective masses. The backwards propagating state is 



the N* (J p 



10 1 



io-° 

10" 7 




I modes 
6000 modes 
5000 modes 
3000 modes 
1000 modes 
500 modes 
30 modes 









— all modes 

' — 6000 modes 
- 5000 modes 

— 3000 modes 

— 1000 modes 
500 modes 

' — 30 modes 





6.4 6.6 

t 



Figure 6: Low mode conttibution to the nucleon two-point function on a 4 x 8 lattice. On the right hand 
side the t = 6,7 data are shown without the logarithmic scale. 

In fig. |5] we show the low mode saturation of the nucleon two-point function and the LMA 
improvement of the corresponding effective mass. The low mode contribution for the ,/V* turns 
out to be negative and is not visible on the logarithmic scale. To gain a better understanding of 
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this behaviour we investigated a small 4 3 x 8 lattice volume where we calculated all the 6144 
eigenvectors of the Wilson operator using the LAPACK library. In fig. ^ we display the nucleon 
two-point function Ci ow (f) on such a single configuration for different numbers of eigenmodes. 
For the positive parity state that propagates from t = into the forward direction the low modes 
are dominant and quickly saturate the correct correlation function. Conversely, for the backward 
propagating states even the sign is wrong until over 85 % of the modes are summed up. Eventually, 
after all modes are included, the correct result is obtained. We also tested non-hermitian LMA for 
baryons. This did not solve the sign problem for negative parity states and, as in the mesonic case, 
no improvement over the conventional point-to-all method was found. 
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Figure 7: Effective masses for the A + (J f 
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, left) and A (J f 



, right). 



For the A + and the A a behaviour similar to that of the 
nucleon was found. In fig. ^ we display the effective mass 
plots. The negative parity partners are not well approximated 
by the lowest modes that exhibit the wrong sign. However, 
LMA reduces the statistical errors for the forward propagating 
states. The errors on the fitted positive parity baryon masses are 
reduced by factors of roughly two, using the hermitian LMA 
method, see table ^[ 





m 


«*LMA 


p 


0.604(16) 


0.595(7) 


A+ 


0.717(20) 


0.713(12) 


A 


0.586(17) 


0.592(9) 



Table 2: Fitted baryon masses in 
lattice units for m u = m c [ = m s . 



4. Conclusions 



We confirm that the hermitian low mode averaging technique is an efficient method to reduce 
the noise of two-point functions and of the fitted hadron masses. It works well for negative parity 
mesons and very well for positive parity baryons but it fails completely for the opposite parity 
cases where low mode saturation does not set in at all. Non-hermitian low mode averaging was 
inefficient for all particles studied, due to its slower saturation behaviour. 

It should be pointed out that computing the low lying eigenmodes is computationally demand- 
ing. Preliminary tests on a 24 3 x 48 volume with m K rj 420 MeV (m n L « 4.8) show that it is 
more efficient to increase the number of point sources than to compute the eigenvectors, to reduce 
the errors on meson masses. For baryons the LMA is cost-efficient. If one further decreases the 
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quark mass then on one hand the costs of computing propagators will drastically increase while 
on the other hand low mode saturation may set in faster, making this method even more efficient. 
Moreover, eigenvectors may be recycled for deflating the Dirac operator, in particular if multiple 
source points per configuration are used, see e.g. refs. |0 In this case the overhead of LMA is 
negligible. 
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